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Abstract. The resonant cyclotron absorption model (Hoshino et al 1992, Gallant and Arons 1994) is very successful 
in describing particle acceleration in plerions. A sensible prediction of this model is the presence of a substantial 
amount of relativistic protons in pulsar winds. Although difficult to detect, these protons may show up through 
their interactions either with the photons in the plerion environment or with the thermal gas in the supernova 
ejecta. Inelastic proton-proton (p-p) collisions are expected to be very effective in young objects, resulting in a 
copious production of neutral and charged pions. Charged pions produced during the first few hundred years after 
the supernova explosion may have time to decay into muons, whose subsequent decay may provide an additional 
source of electrons and positrons in these nebulae, that sums up to the pulsar input. These secondary leptons 
evolve just as the pairs in the pulsar wind, and signatures of their presence could be found, in principle, even in 
the synchrotron spectrum of older objects, pp collisions may remain fairly efficient even in moderately old objects 
resulting in the production of TeV 7-rays and neutrinos. We apply our calculations to the case of the Crab Nebula, 
the best studied plerion insofar, and find that existing data already allow to infer interesting constraints on the 
physical properties of the Crab pulsar wind. 



1. Introduction 

The presence in pulsar winds of a baryonic component co- 
existing with the leptonic one is still an open question. 
According to the most successful model so far for par- 
ticle acceleration at the pulsar wind termination shock 
(Hoshino et al 1992 Gallant and Arons 1994), baryons 
should not only be present in the wind, but even carry 
most of its energy. This makes it mandatory to investigate 
all possible ways in which signatures of their presence can 
be found. 

Plerions result from the interaction of a pulsar with its 
surrounding medium. According to the standard picture, 
pulsars lose most of their energy in the form of a magne- 
tized relativistic wind, whose material component is made 
mostly of electron-positron pairs. The wind starts out cold 
and with a highly relativistic expansion velocity. At some 
distance from the pulsar the wind then needs to be slowed 
down in order to match the boundary conditions imposed 
by the surrounding supernova ejecta expanding at subrela- 
tivistic speed. This happens at a termination shock where 
most of the outflow bulk energy is randomized. Here is 
where particle acceleration presumably occurs, giving rise 
to the non-thermal particle spectrum that is the source of 
the plerion synchrotron and Inverse Compton radiation. 

According to both pulsar magnetospheric theories and 
the available data, the termination shock is most likely 
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a relativistic transverse shock, meaning that the magnetic 
field is almost perpendicular to the flow velocity. Standard 
shock acceleration mechanisms (i.e. diffusive shock accel- 
eration or Fermi I process, and shock drift acceleration) 
are generally believed to be ineffective at this kind of 
shocks (e.g. Hoshino et al 1992). 

An acceleration mechanism that is able to overcome 
the difficulties of the former two is based on resonant 
cyclotron absorption (RCA) in a baryon-loaded plasma. 
Here we briefly review it in order to make clear the as- 
sumptions we make on the proton spectrum and its rela- 
tion with the basic wind parameters. 

Within the framework of the RCA model (Hoshino et 
a/ 119921 Gallant & Arons I1994f) , the ions present in the 
wind have a Larmor radius r Lji which is Am p /Zm e times 
larger than that of the pairs, being both species cold in 
the wind rest frame. When the heavier particles meet the 
shock in the pairs, this appears to them as an infinitesi- 
mally thin discontinuity in the magnetic field. The sudden 
enhancement of the magnetic field is expected to cause 
them to start a bunched gyration, accompanied by the 
collective emission of cyclotron waves. The fundamental 
frequency of these waves is the ion Larmor frequency, but 
also higher harmonics are excited up to the pair gyration 
frequency. The latter can be resonantly absorbed by the 
pairs and cause their acceleration, while the ions slowly 
thcrmalize after crossing the shock. Efficient acceleration 
requires the ions to be energetically dominant in the wind, 
although the pairs are dominant by number. 



2 



Amato et al.: Signatures of high energy protons in pulsar winds 



After crossing the shock the ions drift towards the 
outer parts of the nebula with a velocity of order cr^JR^ 
and slowly thermalize to a temperature that is equivalent 
to a Lorentz factor T of order of the initial wind Lorentz 
factor. The full thermalization takes several gyration pe- 
riods, so, depending on the initial wind Lorentz factor it 
may never be achieved within the nebula. 

While drifting towards the edge of the nebula the ions 
may interact with the plerion radiation field and with the 
thermal material of the ejecta, which partly penetrates the 
synchrotron bubble due to Raylcigh- Taylor instabilities. 
Both these interactions would lead to pion production. 
The products of the subsequent pion decays could lead 
to observable signatures if the interaction rate were high 
enough. The most direct signature would be the instanta- 
neous emission of high energy photons and neutrinos. 

The latter are potentially the preferential channel to 
look at, since their detection would be unequivocally asso- 
ciated with the presence of protons. In fact, a few recent 
papers have considered the question of neutrino produc- 
tion from protons accelerated in pulsar magnetospheres. 
Bednarek & Protheroe l|1997[l computed the flux of high 
energy neutrinos and 7-rays that would be expected from 
nuclear collisions of protons deriving from the decay of 
neutrons produced by the photodisintegration of acceler- 
ated nuclei in a pulsar magnctosphere. Beall & Bednarek 
( 2002 ) computed the neutrino flux that would be produced 
by high energy protons accelerated in the magnetosphere 
of a young fast rotating pulsar surrounded by a super- 
nova envelope, and Bednarek ( 2001 ) considered the effect 
of this contribution on the overall neutrino background. 

None of these previous investigations explicitly consid- 
ered the protons that should be part of the pulsar wind if 
the RCA model for shock acceleration is actually at work 
in plerions. These are exactly the protons we consider in 
this paper, with the aim of finding constraints from obser- 
vations on the characteristics of the pulsar wind. In ad- 
dition to computing the multi-TeV neutrinos and photon 
fluxes we also worry about the secondary electrons and 
positrons produced by ir^ decay. If protons are present 
in fact, this process acts as an extra source of pairs in 
the plerion which is expected to be especially active when 
the plerion is young and compact, but could in principle 
leave relatively long-lasting traces in the nebular radiation 
spectrum. 

In the following we study the efficiency of the pion pro- 
duction processes as a function of time during the plerion 
evolution. We approximate the plerion as spherical and 
expanding at a constant velocity. For the evolution of the 
nebular non-thermal emission we use the Pacini & Salvati 
l|1978J) model (PS hereafter) for a homogeneous plerion. 

2. Pion production and energy losses 

Charged and neutral pions can be generated either by 
photo-pion production off the photons in the nebula or 
by inelastic p-p scattering against the thermal gas in the 
ejecta. While the decay of neutral pions results in 7-ray 



production, the decay of charged particles results in neu- 
trinos and electron-positron pairs (we will refer to both as 
electrons). 

We describe the efficiency of both photomeson produc- 
tion and p-p collisions through a parameter, /"", defined 
as the probability for a single relativistic proton to create 
a pion before losing energy through other kinds of inter- 
actions: 



f = min 



(1) 



where iufc is the proton life-time against losses different 
from p-7 and p-p, and V is the timescale for the process 
of pion production. The luminosity output in pions will 
be: 



(2) 



with L p the wind luminosity in protons. 

In the presence of an intense magnetic or radiation 
field, charged pions may lose a considerable fraction of 
their energy through synchrotron or Inverse Compton 
emission before decaying into muons. The same could hap- 
pen to muons before decaying into electrons and neutri- 
nos. If any of these effects were relevant, this would lead to 
a suppression of the neutrino and electron fluxes at high 
energies. We shall define the quantities / M and f v as: 
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where tj ccay = 2.6 x lO" 8 ^ s and ^ ccay = 2.2 x lO" 6 7/1 s. 
Since the pion energy is roughly evenly shared between 
the final products of its decay, we can write the number 
of neutrinos, electrons and 7-rays injected in the nebula 
per unit time and unit energy interval as: 



J 7 (e 7 ) 
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(4) 
(5) 
(6) 



An important point to be noticed in the above set of equa- 
tions is that in writing Eq.Q) we take into account only 
muon neutrinos, since these are the particles that neutrino 
telescopes are able to detect. The possible effects of neu- 
trino oscillations on the results of our calculations will be 
discussed in section l3~31 

2.1. Photo-meson production 

In the case of photo-meson production the target for high 
energy protons is the plerion emission. The fractional en- 
ergy loss rate of a proton with energy E p (= T m p c 2 ) due 
to pion production results in (Waxman & Bahcall, 1997): 
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2.5 x 10 
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^kpc 



F^mJy] yr- 1 . (7) 



Here we have treated the plerion as homogeneous and 
used the fact that the photon spectrum is a power law, 
F u oc v~ v . We have also made the approximation that the 
main contribution to pion production comes from photon 
energies e 7 « e pC ak=0.3 GeV, where the p-7 cross sec- 
tion peaks due to the A resonance. The numerical values 
are obtained using: cr pC ak = 5 x 10~ 28 cm 2 , £ pca k =0.2, 



Ae=0.2 GeV, v D 



e P eak/(r h) and f3 p ~ 1. 



2.2. Inelastic nuclear collisions 

The energy loss-rate of a relativistic proton due to inelastic 
nuclear collisions can be estimated as 



(8) 



where n t is the target density, <r = 5 x 10 _26 cm 2 and 
( ~ 20% is the average fraction of energy lost by the 
proton. 

As we mentioned in the introduction, the estimate of 
the effective target density, n t , is the main uncertainty in 
our calculation. The distribution of the thermal material, 
in all cases when it is observed, is actually far from uni- 
form. The matter of the ejecta penetrates the synchrotron 
nebulae due to Ray leigh- Taylor instabilities, giving rise 
to concentrations in the form of "filaments". The total 
amount of material contained in these filaments is not easy 
to estimate with confidence. In fact the possibility that the 
gas detected through thermal emission hides some colder 
and denser concentrations cannot be excluded. A compa- 
rably difficult task is to estimate, from the projection on 
the plane of the sky, the fraction of the nebular volume 
that these filaments occupy, their so called "filling factor" . 

Moreover the effective target density for relativistic 
protons, which is our reason for worrying about the dis- 
tribution of thermal matter in plerions, is also affected by 
the intensity and structure of the magnetic field in and 
around the "filaments" , which are not known. The mag- 
netic field affects in fact the physics of proton propagation 
and could in principle be an obstacle to the penetration 
of protons in the thermal gas, leading to a decrease of the 
value of nt, or conversely be the cause of proton trapping 
in the filaments, resulting in the opposite effect. 

We shall summarize all these unknowns in a single pa- 
rameter fi, which we also take to be time-constant. We 
define the effective density of material in terms of the av- 
erage gas density in the plerion as: 



n t — (ith = A* 



m 



= 10 -53- cm 



(9) 



where is the mass of the ejecta contained in the fila- 
ments, which we also take to be time-constant. Using Eq. 
© in Eq. JHJ| we obtain: 



pp 



7 M N0 _ x 
10 ft -^3- yr 



(10) 



which, when compared with Eq. J7J) indicates nuclear col- 
lisions as the main mechanism for pion production in ple- 
rions, if /i is not very much less than 1. 

2.3. Energy losses 

The energy distribution of particles after injection is 
determined, in a plerion, mainly by the following loss- 
mechanisms: escape and expansion losses, and radiative 
cooling through synchrotron and Inverse Compton emis- 
sion. We approximate the plerion as a homogeneous spher- 
ical bubble filled with a uniform magnetic field and ex- 
panding at constant velocity v, and estimate the time- 
scales of the different processes under these assumptions. 
The rate of adiabatic cooling is simply: 



t~ 



R 



\ ) R P c 



yr 
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where R pc is the nebular radius in parsecs and V3 is the 
plerion expansion velocity in units of 10 3 km/s. 

The rate of escape of a particle of mass m s and Lorentz 
factor r s , whose Larmor radius in the nebular magnetic 
field is Vl, can be estimated, based on VB drift, as: 



= ^-3x lO- 3 



P. 2 



B_4 R 2 . m, 



y r 



(12) 



where B-4 is the nebular magnetic field in units of 10 -4 
G and is the particle Lorentz factor in units of 10 6 . 
The rate of synchrotron cooling is: 



*sync = 3 °T C 



7 x 10~ 14 T s 6 B 2 _ 4 
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U B 



yr 



(13) 



Finally, the cooling rate due to Inverse Compton 
Scattering (ICS hereafter) is 



ICS 



-arc 



S u dv 



m s c 
9 , 



rnb 



32 



-v 



KN 



Rat 
c 

Si/ 



dv 



(14) 



where ^kn = m s c 2 /V s /h is the transition frequency be- 
tween interaction in the Thompson and Klcin-Nishima 
regimes, C/ cm b is the energy density of the cosmic 
microwave-background photons, and S v is the volume 
emissivity of the plerion. 

3. Proton signatures during the plerion history 

3.1. Plerion evolution 

All the rates in the previous section depend on quantities 
that change with time during the plerion evolution. We 
determine the dependency on time of the different quan- 
tities within the framework of the PS model. We consider 
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the spin-down power of the pulsar to be continuously con- 
verted into injection of fresh particles and magnetic field 
energy. The spin-down luminosity of the pulsar can be 
parametrized as: 



L(t) 



(l + t/r)« 



(15) 



where the characteristic spin-down time r is a few hundred 
years, and the index a is expected to be around 2 if the 
breaking is purely electromagnetic. 

Treating the plerion as homogeneous and expanding 
at constant velocity, the magnetic field evolution, result- 
ing from the competition between the continuous pulsar 
supply and adiabatic losses, can be determined analyti- 
cally. For times shorter or comparable with r, one can 
make the approximation: 



B(t) 



3Lb 

vH 2 



(16) 



where Lb is the fraction of the pulsar luminosity that goes 
into magnetic energy. 

The other quantity we need to know as a function of 
time in order to estimate ICS losses and the rate of photo- 
meson production is the plerion volume emissivity S v , 
due to non-thermal radiation. To this purpose, we need 
to determine the time-dependences of the spectral den- 
sity N(E,t) of the emitting particles (pairs). For the wind 
particles we assume that they are injected with an initial 
spectrum in the form of a broken power-law: 



ME) 



K { {E/E\) 



E mm < E < E\ 

Ei < E < E max 



(17) 



where K w is related to the fraction of pulsar luminosity 



that is converted into pairs, L pa irs through 



K w {t) 



L 



pairs 
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Ei 



2-7i 



72 



'EL 



Ei 



(18) 



In the following we shall assume E\ to be constant with 
time and independent on the wind Lorentz factor: we take 
Ei = 1.5 TcV, roughly corresponding to the energy of 
the particles radiating at the frequency at which the high 
energy break of the spectrum of the Crab Nebula is ob- 
served. As to -E max , this will depend on both time and 
the wind Lorentz factor, due to the dependence on the 
latter of the proton energy and on the dependence on the 
former of the nebular magnetic field strength. In fact we 
assume the maximum energy to which the pairs can be ac- 
celerated to be given by E max = mm(E p , E cut ), with E cut 
determined by the condition T acc (E cut ) = T lossos (E cut ), 
i. e. the acceleration time cannot exceed the time-scale for 
energy losses. Assuming T acc rj fl^ 1 (with O the particle 
Larmor frequency) and the losses as due to synchrotron 
emission in the average nebular magnetic field, one finds 
E cut {t) = m e c 2 ^&ne/{a T B{t)). 



We mentioned that charged pion decay provides an 
additional source of pairs in the nebula. Their evolution is 
completely analogous to that of the wind particles, so that 
apart from the different initial spectrum, the treatment 
below applies to both components. 

The spectral distribution at a time t, N(E, t), is found 
as a function of the injection spectrum J = J w + J e (see 
Eq. |J5J), as: 



N(E, t) 



dt 

J{Ek,U) g^(E,Ei,t) dE t 



(19) 



where U is the time at which the particle was born and Ei 
its initial energy. N(E, t) is readily computed once dU/dE 
is known as a function of E and t. 

The needed relation between E, Ei and U may be 
found from the equation describing the single particle en- 
ergy evolution due to synchrotron and adiabatic losses: 
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where ci = (1/6tt) <jt c/(m e c 2 ) 2 is the constant relat- 
ing the synchrotron power to magnetic field strength and 
particle energy. For dU/dE we find: 



dU 
dE 
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2 E,, 
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(21) 



where we have used the definition of the so called break 
energy: 

2 



E b (t) 



a t B 2 (t) 



(22) 



For times shorter than r, when B(t) oc t 1 (Eq. (|16fl ). 
so that Eb(t) oc t, one finds: 
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2(1 + E/E b (t)) E t 
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(23) 



and the solution of Eq. (|19|) can be approximated as: 



N(E,t) 



2 t 



J(E,t ), , E<E b {t) 



7i + 72 I (E b (t)/E)J(E,t ) , E b (t) < E 



(24) 



with to ~ for t < r and t ~ r afterwards. 

The evolution of the wind particle spectrum as a func- 
tion of time is plotted in Fig. for the case T = 10 6 , 
7i = 1.5, 72 = 2.2 and L = 3 x 10 39 erg s" 1 , with 
Lb = ipairs = L p = (1/3)L. We plot the full solu- 
tion of the problem taking into account the real time- 
dependences arising from the varying pulsar input (Eq. 
(|15fl with r = 700 yr) together with the approximate ex- 
pression in Eq. I|24|l . 

Once N(E, t) and B{t) are known, the volume emissiv- 
ity S v that enters Eq. (|14(l is computed from synchrotron 
theory as: 



ciB 



2c 2 4 7T R% V c 2 B 
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(25) 
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A t=10 yr «t=500 yr 



30 



E[eV] 

Fig. 1. The time evolution of the wind particle spectrum in 



the case 71 = 1.5, 72 = 2.2, L p 



10 



39 erg s 1 



The 



different curves refer to different epochs after the supernova ex- 
plosion: t = lyr (solid), t — lOyr (dotted), t — lOOyr (dashed), 
t — 500yr (dot-dashed). The points trace, at each time, the 
approximation given in the text. 



where C2 = 0.3 e c/(2 ir (m e c 2 ) 3 ). Knowing Sf, ync one 
can also compute the ICS emissivity of the plerion, Sl cs , 
following, for instance, Blumenthal & Gould 1|197(J[) . From 
the total emissivity S v = Sl ync + Sl cs we also obtain the 
flux F v that appears in Eq. Q: 



F v (y) 



3% *M 



(26) 



3.2. Pion production efficiency 



The time-evolution of the magnetic and radiation field in 
the nebula as derived in the previous section was used to 
compute the time-history of pion production efficiency. 

In the plots below we show the behaviour of / M , 
fZy and fp p as a function of time after the supernova ex- 
plosion for various values of the pulsar luminosity L, of 
the wind Lorentz factor T and of the remnant expansion 
velocity v. In all the plots in this section we assumed the 
pulsar energy input to be evenly shared between magnetic 
field, protons and pairs: Lb = L p = L pa j rs = (1/3) Lq. 
To compute the radiation field, we assumed the injection 
spectrum of the pairs as in Eq. (|17(l with 71 = 1.5 and 
72 = 2.2. We consider a monoencrgctic distribution of 
protons corresponding to the Lorentz factor of the wind. 

In Fig. J21 we plot / A '(£ , 7r ,i) for the value of E v cor- 
responding to the highest energy pions produced for each 
of the values of T considered. 

It is apparent from the figure that for Lorentz factors 
up to 10 5 pions never suffer any severe energy losses, not 
even at early times, when the magnetic and radiation field 
in the plerion are most intense. For energetic enough pi- 
ons, instead, radiation losses at early times do become 
important, decreasing the decay probability by a factor 
that can be of order 10 2 for the lowest considered value 
of the expansion velocity and highest considered pulsar 
luminosity (filled triangles). 




Fig. 2. The time evolution of / M for different values of the 
pulsar luminosity, of the remnant expansion velocity and of the 
wind Lorentz factor. The curves plotted in different panels are 
computed for different values of the wind Lorentz factor. The 
different curves in each panel are marked as explained below. 
Empty triangles: L[erg/s] = 3 x 10 39 , «[km/s] = 10 3 . Empty 
squares: L[erg/s] = 3 x 10 39 ,v[km/s] = 10 4 . Filled triangles: 
L[erg/s] = 3 x 10 40 ,«[km/s] = 10 3 . Filled squares: L[erg/s] = 
3 x 10 40 ,u[km/s] = 10 4 . The vertical lines represent the time at 
which the optical depth for neutrinos drops below 1 in the case 
when v = 10 3 km/s (solid) and when v = 10 4 km/s (dashed). 




Fig. 3. The time evolution of /" for the highest energy pions 
corresponding to different values of the pulsar luminosity, of 
the remnant expansion velocity and of the wind Lorentz factor. 
Same notation as in Fig. 



Losses at early times are found to be due mainly to syn- 
chrotron emission and their dependence on v and L arises 
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from the relation between the magnetic field intensity and 
these two quantities (Eq. l(TBJ). It is interesting to no- 
tice that in all the cases considered, / M becomes 1 at very 
early times, and actually earlier than the time at which the 
remnant stops being optically thick to neutrinos, which we 
mark as a reference time in Fig. This reference time 
is computed considering the fact that the main source of 
optical depth for high energy neutrinos are nuclear colli- 
sions that give rise to muons: v + X — > fx + Y, The cross 
section for this process (e.g. Learned & Mannheim 2000) 
is a function of the neutrino energy alone at the energies 
we consider. Hence, the time at which t v = Rn n t o v x 
drops below unity for the highest energy neutrinos cor- 
responding to each value of the Lorentz factor we con- 
sider depends only on the Lorentz factor itself and on the 
density evolution of the remnant (i.e. on its expansion 
velocity, once the total mass is fixed and assumed as uni- 
formily distributed). One has: r„ oc T°/v, with f3 = 0.5 if 
F < 4 x 10 5 and (3 = 0.2 otherwise. 

The same legenda as for Fig. J2J applies to Fig. ©, 
where the evolution of f v is shown (again in the case when 
it is minimum, i.e. for the maximum muon energies). 

As we could expect, radiation losses for muons are 
more severe than for pions, and for the most energetic we 
are considering (lower right panel) j v actually becomes 1 
later than / M . 

It is important to notice that as far as only the non- 
thermal radiation field is taken into account, synchrotron 
losses are always dominant with respect to losses due to 
ICS, according to our calculation. The latter can become 
important, though, if also the thermal radiation field is 
taken into account, but this effect should be limited to 
earlier times than that at which t v drops below 1 (Beall 

& BednarekEOnil)- 

As far as f* is concerned, we computed the efficiency of 
both photo-meson production and p-p collisions. For the 
latter process we considered 1 Mq of material uniformily 
distributed within the remnant and used \i = 1 in the 
calculations. The results are shown in Fig. (J3J. 

We find that radiation losses are never effective to limit 
the efficiency of p-p collisions, not even in the very intense 
initial magnetic field of a bright, slowly expanding plerion. 
Conversely, for high enough values of the wind Lorentz 
factor r and low enough values of the wind expansion 
velocity, synchrotron losses are the main limit to photo- 
meson production. However, we can see from Fig. El that 
photomeson production never reaches efficiencies higher 
than about 30%, at least for the prototypical spectrum 
that has been assumed here. 

It should be noticed, finally, that in the plots shown 
above for / w , / M and f u , when computing both ICS losses 
and photo-meson production we have only included the 
photon spectrum deriving from synchrotron and ICS emis- 
sion of the wind particles. In doing that we have neglected 
the thermal emission from the gas and dust in the fil- 
aments and the radiation deriving from pion decay. The 
consistency of the latter approximation was verified a pos- 
teriori: after calculating the radiation flux from the prod- 




Fig. 4. The time evolution of f for different values of the 
pulsar luminosity, of the remnant expansion velocity and of the 
wind Lorentz factor. The solid curves refer to the efficiency of 
nuclear collisions while the dashed ones show the efficiency of 
photo-meson production. Curves referring to different values 
of the pulsar luminosity and remnant expansion velocity are 
marked according to the same notation as for Figs. (5J and 



ucts of pion decay we checked that adding this contribu- 
tion to the target radiation for ICS losses did not change 
the results plotted in Figs. 

3.3. Photon and neutrino fluxes from a "young Crab" 

In this section we use the theory introduced above to com- 
pute both the electromagnetic and neutrino emission from 
a young plerion with Crab-like parameters, as a func- 
tion of time after the supernova explosion. We consider 
a pulsar whose initial spin-down power is L — 3 x I0 39 
erg/s. We assume, as in section EOl that L is evenly con- 
verted into electron, proton and magnetic field energy: 
Lb = L p = Lpairs = (l/3)Lo. The distance of the pul- 
sar and the plerion expansion velocity are taken to be the 
same as for the Crab, and 1 Mq of thermal material uni- 
formily distributed in the plerion volume is considered as 
the target for p-p collisions. 

The time evolution of the pion production efficiency in 
such a plerion was computed in section l3~^l (curves marked 
with empty squares in Figs. @-(QJ). Taking the results 
from that section we estimate the global photon spectrum 
of the plerion. The high energy part is entirely due to tt° 
decay, while at low energies we estimate the synchrotron 
emission of both wind and secondary pairs together with 
the ICS emission of the former. The overall spectrum at 
different times during the plerion evolution is reported in 



Fig. ©. 
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It is apparent from the figure that, for the set of pa- 
rameters used, the radiation due to the secondary pairs 
could only be detectable over the emission due to the wind 
particles, at early times (up to few tens of years after the 
supernova explosion) and at photon energies in the range 
100 MeV-1 TeV. 

We notice that 7-rays from ir° decay should be de- 
tected even by present TeV detectors during the first few 
tens of years after the supernova explosion, if a similar ple- 
rion was born within the galaxy. The ir° decay 7-ray flux 
always exceeds the ICS emission at high enough energies 
for times upto several tens of years at least. 
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Fig. 5. Electromagnetic spectrum of a "young Crab" at dif- 
ferent times after the supernova explosion. The thick curves 
correspond to the synchrotron (at low energies) and ICS (at 
higher energies) emission of the wind pairs, while the lighter 
curves refer to the emission of the pion decay products. The dif- 
ferent panels correspond to different values of the wind Lorentz 
factor r and the different line types refer to different times. 



In Fig. 10 we show the expected neutrino flux from 
the same plerion. We hnd that this is generally well above 
the atmospheric background of a detector with an angular 
resolution of 2° x 2°. Moreover a galactic object described 
by our set of parameters could give a flux of neutrinos 
already detectable by telescopes like AMANDA II for the 
first ten years. It is clear from the figure that the proba- 
bility of detecting such objects will increase considerably 
with the upcoming km 2 neutrino telescopes, since their 
lower threshold allows neutrino detection even from older 
plerions, if the pulsar wind contains a non-negligible frac- 
tion of protons. 
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Fig. 6. Neutrino spectrum of a "young Crab". The notation 
is the same as for Fig. ©. The solid thick line is the detection 
threshold for ICE-cube CHill l2UUll . while the dot-dashed line 
is the atmospheric neutrino background in a 2° x 2° bin. 



4. Constraining the parameters of the Crab pulsar 
wind 

In this section, we are going to use the Crab Nebula, the 
prototype plerion, to discuss how the information coming 
from high energy photon and neutrino detection could be 
used to constrain the physics of plerions. 

The parameters playing a role in determining the flux 
of high energy photons and neutrinos are all quite well 
known for the Crab Nebula, except for the parameters that 
enter in the determination of /i, the wind Lorentz factor 
r, and the fraction of the pulsar spin-down power that 
goes into acceleration of protons. The relevant parameters 
for our purposes are, from the equations in section [21 the 
plerion age t = 950yr, its distance d ~ 2 kpc, its expansion 
velocity, V exp — 1.5 x 10 3 km s _1 , the average magnetic 
field intensity B ~ 3 x 10~ 4 G, and finally the amount of 
mass contained in the filamentary shell. For the latter, the 
most recent estimates give Mm = (4.6 ± 1.8)Af Q (Fesen 
et a^ ll997|) . We shall consider a value Mm = 1M Q and 
values of // in the range 1 < \i < 20, according to the 
upper limit on n t derived by Atoyan & Aharonian ( 1996), 
based on the compatibility between the bremsstrahlung 
7-rays and the observed flux from the Crab Nebula above 
1 GeV. We notice that if the mass estimated by Fesen et 
aL( 1992J, were uniformily distributed in the nebula, this 
would correspond to /i = 5. 

Inserting the numbers just quoted into the equations in 
sectionFJJwe can estimate all the relevant time-scales. One 
can readily see that what is important for computing f n is 
just the comparison between t osc and i oxp (whose relative 
importance depends on T) and t pp , since the other three 
time-scales are much longer than these. From Eqs. (TJJ, 
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<jTU|) . l(TT|l and (U2J, we have 
f = min 



= min 



-, ^cxp . / -, ^esc 

1, - — mm 1, 

£pp \ ^exp , 

1,2.8 x 1CT 5 // min (1, — 



(27) 



where Tq is the wind Lorentz factor in units of 10 6 . 

For the same values of the parameters the muon losses 
due to synchroton emission or ICS (Eqs. lflH|) and iQ ) 
are completely negligible on the time-scale of their decay, 
even for the highest energy we consider. The same is true 
for pions, whose proper decay times are about a hundred 
times shorter. 

In Figs. Q) and ©, we plot, respectively, the expected 
neutrino and photon fluxes from the Crab Nebula as a 
function of energy, for different assumptions on the wind 
Lorentz factor and different values of the parameter /i. 
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Fig. 7. The flux of neutrinos expected from the Crab Nebula 
as a function of energy, for different values of the parameter 
/x and for different values of the wind Lorentz factor F. The 
different curves in each panel correspond to fj. = 1 (solid), 
H = 5 (dotted), fj, = 10 (dashed) and /i = 20 (long-dashed). 
The considered values of T range from 10 4 to 10 7 , as specified 
in each panel. The solid thick line is the detection threshold 
for IceCube (Hill 2001), while the dot-dashed line is the atmo- 
spheric neutrino background in a 2° x 2° bin. We assume 60 
% of the pulsar wind energy to be carried by protons. 



We have assumed that protons are the main energy 
carriers in the wind so that we put in them about 60% 
of the pulsar spin-down luminosity L p — 0.6 Libt — 3 x 



10 erg s , being Lxot = 5 x 10 erg s . Our choice 
of the energy fraction into baryons comes from the fact 
that this is the percentage of the pulsar power which is 
required to be in baryons if the luminosity contrast of the 
Crab Nebula wisps is to be explained through magnetic 



Fig. 8. The flux of high energy photons expected from the 
Crab Nebula, for different values of the parameter n and of 
the wind Lorentz factor V. The notation is the same as in Fig. 
Q. The expected fluxes are compared with the available data 
and upper limits (de Jager et al 1996 , Aharonian et al 2000 
and references therein). Also shown as a thick curve in each 
panel is the ICS flux computed based on this model. 



field compression due to the ions gyration in the shocked 
background pair plasma (Spitkovsky & Arons 1999). 

The neutrino and photon fluxes were computed from 
Eqs. (0} and © respectively, with the pion differential 
spectrum per unit time computed in the scaling approx- 
imation for the p-p scattering cross-section fBlasi l2T)01fl . 
In the scaling regime a monoenergetic proton distribution 
with energy E p leads to the injection of pions described 
as: 



dE^dt 

where the function g„ is 
^(^(l-^ + e- 18 * 



/1.34 , 



and is found from the condition 



f n Lp — 
leading to: 



E^dE^ 



g v {x)dx 



(28) 



(29) 



(30) 



(31) 



It is apparent from Fig. Q that, for the generally 
favoured value of T ~ 10 6 (Kennel & Coroniti I1984al 
and I1984bl Gallant & Arons 1994), and for the value of 
L p ~ 0.6Ltot favoured by Gallant & Arons (1994) and 
Spitkovsky & Arons l|1999|l , the flux of neutrinos from the 
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Crab Nebula would be well above the atmospheric neu- 
trino background already for /i = 1. The chances of de- 
tecting multi-TeV neutrinos from the Crab Nebula with 
IceCube look promising, based on this figure. 

However, part of the section of the parameter space 
that would lead to a detectable flux of neutrinos can al- 
ready be excluded based on the available high energy 7-ray 
data. In Fig.© we plot the flux of high energy photons 
expected from ir° decay. It is clear that in some cases the 
computed fluxes at multi-TeV energies are above the mea- 
surements and upper limits. The latter, on the other hand, 
can be comfortably explained, whithin the error bars, as 
due to ICS (de Jager et »/fTT)9(i Atoyan & Aharonianl.1996 
Aharonian et aZ l2U(J(J|) . Even within the framework of our 
crude model, which not only treats the plerion as homo- 
geneous but also neglects the contribution to the target 
photon field of the ambient thermal emission, the com- 
puted ICS flux is only slightly below the TeV data points. 

In the following we shall review our expectation for the 
flux of neutrinos from the Crab Nebula in the light of the 
constraints arising from 7-ray observations. 

Before discussing these limits, however, we should also 
consider the other possibility we mentioned to put con- 
straints on the three unknowns /x, L p and T. This is of- 
fered, at least in principle, by the signatures that may be 
left by the secondary pairs in the synchrotron radiation 
spectrum of the nebula. These pairs would be injected in 
the Nebula according to Eq. (J5J, but the injection rate (de- 
scribed by Eq. I|28|) and depending on /*" and L p through 
Eq. ABU ) has to be considered now as a function of time. 
Since we are mainly interested in low energy pairs, with 
long synchrotron life-times, we need to take into account 
the whole history of pion production in the nebula. This 
is done as discussed in sections l3.ll and [3.21 

We show in Fig. the comparison between the Crab 
Nebula observed radiation spectrum and the synchrotron 
flux from secondary pairs computed in the case when 
L p (t) = 0.6 L Tot (t), with L Tot (t) described by Eq. lfT5|) . 
where r = 730 yr and a — 2.3 (Groth 1975). In the same 
figure, we also show the contribution arising from 7r° de- 
cay, in the 7-ray part of the spectrum. It is apparent that 
the emission produced by the secondary pairs, according 
to our model, is a completely negligible fraction of the 
present day synchrotron emission of the Crab Nebula and 
cannot provide any additional constraint on the parame- 
ters we are interested in. 

In Fig. I|10fl we use the information that can be ex- 
tracted from high energy data (Fig. JSJ) to restrict the 
allowed parameter space for the values of Lp/Lxot , T and 
/1. The allowed regions, for each choice of /x are those lying 
below the corresponding curve, as explained in the figure 
caption. 

It is interesting to notice that, if // = 10 or larger, the 
possibility that 60% of the pulsar wind energy might be in 
protons with a Lorentz factor T — 10 6 would be excluded 
by the data. 

We now use the constraints in Fig. (|10f> to estimate 
the maximum number of neutrino events that one could 
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Fig. 9. Consequences of pion production on the electromag- 
netic spectrum of the Crab Nebula. The thick curve and the 
points represent the observed emission, while the thin curves 
are fluxes from our model. At low frequencies the model flux is 
due to synchrotron emission associated to secondary electrons 
and positrons, while the contribution at high frequencies is the 
same shown in Fig. © and comes from 7r° decay. The different 
panels refer to the emission computed for different values of 
the wind Lorentz factor while the different line-types are asso- 
ciated to different values of the parameter fj, as specified in the 
each panel. 
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Fig. 10. The allowed parameter space for T and L p / L^ot as 
constrained through the comparison with observations of the 
photon fluxes from pion decay computed within the framework 
of our model. The allowed regions for each of the considered 
value of jj, are the ones lying under the curve of corresponding 
type: fi = 5 — » solid, fj, = 10 — > dashed, /1 = 20 — > dot-dashed. 
No constraints can be put on the wind parameters for n = 1. 

expect to detect from the Crab Nebula in 1 yr of observa- 
tion with a detector of area 1 km 2 , such as IceCube will 
be. The result of this calculation is shown in Fig. I|llf> 
where we plot the maximum number of muons per year 
as a function of the wind Lorentz factor. N p is computed 
taking into account the probability for a neutrino to inter- 
act with the matter in the detector and create a muon. For 
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each value of T we consider the maximum allowed value of 
f n L p with the requirement that /i < 20 and L p < 0.9ixot- 
This determines J v {e u ). iV„ is then found as: 



lyr x 1km 
4 7T d 2 



(32) 



ITcV 



where the maximum neutrino energy is eMax = m p c 2 T/4 



6 J 



and P vp , = 1.3 x 10""e 



.,,. v with 13 = 1 if 1 TcV < e v < 
100 TeV, and (3 = 0.5 otherwise (Gaisser et al\]MB- 
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Fig. 11. The maximum number of muons per year that would 
be produced by neutrinos coming from the Crab Nebula in a 
km 2 detector. Also shown, as a dot-dashed line, are the back- 
ground counts, assuming an angular resolution of (1°) , such 
as IceCube should have. 

The background counts we show in Fig. I|ll[) are 
those computed for IceCube by Alvarez-Muhiz and Halzen 
( 2001 ) at the position of the Crab Nebula. The maximum 
number of neutrino events allowed by the 7-ray data is 
above the background for basically the entire range of val- 
ues of r we considered. The effect of neutrino oscillations 
would be to reduce these fluxes. For the case of maximal 
mixing, equal abundancies of neutrinos of the three flavors 
are obtained, so that the flux of muon neutrinos is about 
half of that calculated above. 



5. Conclusions 

In this paper we have anylized the observational conse- 
quences that the presence in pulsar winds of a proton 
component, energetically dominant, might have. 

We have found that when the plerion is very young 
(up to a few tens of years after the supernova explosion) 
the main loss mechanism for these protons is likely to 
be pion production due to nuclear collisions. A large flux 
of neutrinos, 7-rays and secondary pairs can be expected 
during these early stages. The synchrotron emission from 
secondaries is likely to be completely negligible, while the 
TeV photon and neutrino emission from a Galactic object 
could be detectable even by present- day telescopes during 
the first few tens of years after the supernova explosion. 
Detection of neutrinos from older plerions should await 
the upcoming km 2 detectors. 



We have also estimated the present-day TeV photon 
and neutrino emission from the Crab Nebula. The existing 
TeV observations already exclude part of the parameter 
space, so that not all combinations of T, L p and \x are 
allowed. 

Finally we have computed the maximum number of 
neutrino events allowed by the 7-ray data according to 
our modeling and found it to be larger than the at- 
mospheric neutrino background (as already suggested by 
Arons, 1998). Such a flux is likely to be detected by a 
km 2 detector such as IceCube. Neutrino detection can be 
therefore considered as a promising means to help clinch 
the question of how large a fraction of the pulsar spin- 
down energy goes into relativistic baryons. 
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